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^^ ■ In a mesoscopic metal in proximity with a superconductor, the electronic 

f— ^ . conductance is enhanced in a very energy-sensitive way. In this paper, we 

0\ I discuss the spectral conductance of a proximity superconductor from both the 

^^ ■ theoretical and experimental point of view. The dependence of the spectral 

"jrt ' conductance on the phase-breaking length, gap of the superconductor and in- 

^ I terface transparency is theoretically investigated. We present experimental 

data on the re-entrance of the normal-state conductance at very low tem- 
perature and bias voltage. A complete description of the experimental data 



O . needs taking into account heating of the reservoirs by the bias current. In 



addition, we show that the energy sensitivity of the proximity effect enables 
one to access the energy distribution of the conduction electrons inside a 



K^ , mesoscopic sample. 



1. INTRODUCTION 

1.1. The general context 

In the past few years, the study of the proximity effect in hybrid struc- 
tures made of a normal metal in contact with a superconductor has known 
a remarkable revival .El'Bd This originated from both the availability of a new 
generation of tailored-made submicron-sized samples and the new under- 
standing of the coherence effects in electron transport at the mesoscopic 
scale. The proximity effect is the generic name for the phenomena appear- 
ing at the interface of a normal metal (N) with a superconductor (S). In 
a N-S geometry, the N-metal shows superconducting-like properties includ- 
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ing a magnetic screenings a modification of the density of statea3 and an 
energy-dependent conductance enhancement HQ'Q In a S-N-S geometry, the 
Josephson effect is the manifestation of the coherence in the proximity effect. 
Eventuahy, the conductance of a N-I-S tunnel junction is governed by the 
tunnehng of electron pairs. 

In a N-I-S tunnel junction, the BTK theoryB describes the cross-over in 
the behaviour of the interface conductance as a function of tunnel barrier 
transparency. The interface conductance is twice the normal-state one in 
the high-transparency limit {Gns = '^Gnn) while it vanishes {Gns -^ O) 
in the low-transparency limit. The discovery of the zero-bias anomalyii3 
showed a clear discrepancy with the BTK theory : in a confined geometry the 
conductance of the N-I-S junction does not go to zero at zero temperature 
but increases at low temperature and bias. The physical picture for this 
effect is that the confinement of electrons by the disorder in the vicinity of 
the N-S interface makes the transmission probabilities add coherently. This 
enhances drastically the tunneling of pairs.O'Ej 

Diffusive electron transport in the N metal part of a N-S junction is 
also significantly enhanced by the proximity effect .B'Llcl The observed large 
amplitude of the conductance enhancement used to appear in contradic- 
tion with the prediction of the Bogoliubov-de Gennes equationalj and the 
random matrix theoryJHI The observation of the re-entrance effect has been 
the key experiment that made the con ri e ct i on b etw een the high-temperature 
regime and the low-temperature regime .113113113111 At finite temperature or 
bias voltage, the conductance of the "proximity superconductor" made of 
the N metal in proximity with S is significantly enhanced, while it should 
coincide with the normal-state value at zero temperature and bias. These 
features arise clearly from the quasiclassical theory. 

In this paper, we will focus on the conductance enhancement in very 
small N-S structures. We will consider the regime where the N-S interface 
conductance is large compared to the metallic conductance of the N metal 
wire. The small size of the structures enables us to investigate the trans- 
port properties at temperatures well below the characteristic energy of the 
proximity effect. 

In the following, we will first introduce the main theoretical concepts of 
the proximity effect. We will concentrate on the spectral conductance and 
its dependence on various physical parameters. Afterwards, we will present 
a thorough experimental study of the re-entrance effect in a mesoscopic 
structure made of a normal metal wire in contact with a superconductor. 
We will compare our experimental results with the prediction of the theory 
and conclude about the energy distribution in the reservoirs. Eventually, 
we will show that the large energy-sensitivity of the proximity effect can be 
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used to test the electron energy distribution inside a mesoscopic sample. 

1.2. The Andreev reflection 

Let us consider the interface between a normal metal and a supercon- 
ductor. Well below the superconducting transition temperature of S and 
at low bias voltage, the thermal energy k^T and the electrostatic energy 
eV are much smaller than the gap A of S. In this regime, single electrons 
cannot enter the superconductor because of the absence of avalaible states 
at the same energy. As a consequence, electrons arriving from N will be 
Andreev-reflected at the N-S interface, see Fig.|l]. In this process, the elec- 
tron is retro-reflected as a hole, while a Cooper pair is transmitted in the 
superconductor J12I 

At the Fermi energy, the reflected hole has the opposite spin, the oppo- 
site velocity and the same momentum as the incident electron. The reflected 
hole traces back the trajectory of the incident electron. The reflection of a 
hole may also be seen as the absorption of a second electron. As the Andreev 
reflection correlates the two electron states, it is equivalent to the diffusion 
of pairs in the N metal. We call the diffusing pairs in N " Andreev pairs" as 
their existence is not due to an intrinsic attractive interaction in N, but to 
a remote effect which is the Andreev process at the N-S interface. 

If one looks more precisely, the reversal of every electron velocity com- 
ponent is perfect only at the Fermi levelrl Let us consider an electron with 
a small extra energy e compared to the Fermi level energy Ep. The electron 
wave- vector ke = kp + q 'is larger than the Fermi wave- vector kp. The re- 
flected hole has a slightly different wave- vector k^ = kp — q, see Fig. ||. This 
means that the incident and reflected particles will have a difference 

6k = 2q = kp-^. (1) 

£!jp 

in wave-vector. This difference concerns the wave-vector component which 
is perpendicular to the N-S interface. With an arbitrary angle of incidence, 
the change in wave-vector will change the trajectory angle with respect of 
the N-S interface.E3 In brief, the retro-reflection is imperfect at finite energy. 



1.3. The diffusion of "Andreev pairs" 

In the following, we will focus on the case of a mesoscopic normal metal 
in the dirty limit. This means that the elastic mean free path Ig is much 
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Fig. 1. Left : Schematic of the Andreev reflection process. An incident 
electron with an extra energy e compared to the Fermi energy Ep hits the N- 
S interface from the N side. The reflected hole-like particle and the incident 
electron have a slight wave-vector mismatch 2q. Right : Relevant length 
scales with their schematic respective amplitudes in a metallic thin film. 
The energy-dependent coherence length Lg is the length over which the two 
components of the Andreev pair acquires a phase difference of order vr. 



smaller than the sample length L which is itself smaller than the phase- 
breaking length L^ : Ig <^ L <^ L^,. 

Let us consider the trajectories of an electron incoming on the N-S inter- 
face and of the reflected hole nearly retracing the trajectory of the incident 
electron. At zeroth order, the phase acquired by the electron is eaten up as 
the hole retraces back the trajectory of the incident electron. Looking a little 
closer, the wave-vector mismatch between the electron and the hole makes 
the phase difference between the hole and the electron increase monotoni- 
cally. After diffusion over a distance L from the interface, the phase shift 
between the two particles is of order vr at a distance equal to the energy- 
dependent coherence length : 



L, 




(2) 



D being the diffusion coefficient in N. At the same time, the trajectories 
of the electron and the hole are shifted by a distance of order the Fermi 
wavelengthHj The exact value of the spatial shift depends on the incidence 
angle on the N-S interface. Further diffusion of the two particles will there- 
fore be different and the pair will break apart. In this respect, the length L^ 
actually plays the role of the coherence length of the Andreev pair. 

The coherence length L^ coincides with the thermal lengthEll Lt at 
the energy e = 2TrkBT. The thermal length Lt is relevant as soon as one 
considers the whole electron distribution at thermal equilibrium. This is for 
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instance the case in the Josephson effect. At the Fermi level, the length L^ 
diverges. Indeed, the coherence length of an Andreev pair is limited by the 
phase-breaking length L^p of a single electron. The average coherence length 
of an Andreev pair therefore varies from about the thermal coherence length 
Lt at high energy (e ~ ksT) to the phase-coherence length L^ at low energy 
(e~0). 

The correlation between the electron and the reflected hole can also be 
described in terms of energy with the equivalence : 

Le = L <^ e = ec (3) 

where 

ec = jj (4) 

is the Thouless energy, or correlation energy, of a sample of length L. The 
Thouless energy is the fundamental energy scale for the proximity effect at 
the spatial scale L.Q At a given distance L, only electrons with an energy 
below the Thouless energy are still correlated as Andreev pairs. The other 
Andreev pairs are broken. 



2. THE THEORETICAL BACKGROUND 
2.1. The Usadel equation 

The diffusion of superconductivity in an inhomogenous structure can be 
described by the Gorkov Green's functions. In the limit where all the exper- 
imentally relevant length scales are much larger than the Fermi wavelength, 
the si mplifi c atio n of t his fully quantum theory into the quasiclassical theory 
is valid. E3E30'e3'e3£ZK The Usadel equations are obtained in the diffusive 
regime where the elastic mean free path is small. In this framework, weak 
localization effects and conductance fluctuations are neglected. 

Let us restrict to the case of a one-dimensional N-S structure in zero 
magnetic fleld and no superconducting phase gradient. The latter condition 
is fulfllled as soon as there is only one superconducting island. In the absence 
of electron-electron interaction in N, the Usadel equation in the N metal is : 

hD6le + {2i€--Y cos9}sme = (5) 

The complex proximity angle 9 = 9i + i92 is related to the anomalous 
Green function F : 

F{e,x) = -ism[9{e,x)] (6) 
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The anomalous Green function F is often called the pair amplitude although 
if differs from the actual pair density which takes into account the energy 
distribution of the electrons. 

The complex angle 9 and the pair amplitude F are functions of both the 
energy above the Fermi level e and the distance from the S interface x. The 
Usadel equation features a non-linear diffusion equation for the proximity 
angle 9. Let us note that the phase-breaking length L^ enters the Usadel 
equations as a cut-off for the proximity effect. The inelastic mean free path 
Lin is included in L^, but does not show up by itself. It enters directly in 
the problem only through the energy redistribution of the electrons in N. 

At the contact with an N metal reservoir, the angle 9 is zero. At the 
contact with a superconducting island S, a complete treatment implies solv- 
ing the self-pconsistency equation for the pair potential in the vicinity of the 
interface.E3™ Here we will restrict to the idealized case of a step-like pair 
potential. The pair potential will be considered as being constant and equal 
to the gap in S, and zero in N. With this assumption, the boundary condition 
for a perfectly transparent N-S interface is 9 = iir/l at energy zero or well 
below the gap A. 



2.2. The spectral conductance 

Electron transport is non-local at the mesoscopic scale. In the frame- 
work of the quasiclassical theory, the electron transport is nevertheless de- 
termined by a local conductivity a{e, x) which expresses as : 

cr{e,x) = cJiv cosh^ 92{e,x) (7) 

where ajsf is the normal-state conductivity. From this relation, the local 
conductivity is always larger than the normal-state conductivity. Since the 
function 02 is strongly energy-dependent, so is the conductivity enhance- 
ment. 

To proceed with transport properties, one needs to know the occupa- 
tion of current-carrying states. The non-equilibrium quasiparticle energy 
distribution function can be derived from the out-of-equilibrium Keldysh 
Green function.c3 The generalized quasiparticle distribution function has 
two components : an odd part and an even part. The odd part is related 
to the energy distribution functions of the electrons fe and the holes fh as 
• /o(e) = l/2(/e + fh), and reduces to tanh{e/2kBT) at equilibrium. The 
even part / = l/2(/e — fh) is related to the imbalance in the population of 
holes and electrons, and reduces to zero at equilibrium (no current). It can 
therefore be named the out-of-equilibrium part of the distribution function. 
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At any point in a wire of section S, the spectral current i{e, x) is related 
to the spatial derivative of the distribution function /. In this work, we 
will consider the regime where the inelastic mean free path Lj„ is much 
larger than the sample length : Lj„ ^ L. In this regime, a conduction 
electron keeps its energy while travelling through the sample. Therefore, the 
transport channels at different energies are independent and the spectral 
current z(e) is constant along the wire. This permits to write the spectral 
current as a function of the spectral conductance g{e) and the difference 
A/(e) of the distribution functions / at the two extremities of the considered 
wire : 

i{e)=aie,x)S^^^=gie)Af{e) (8) 

As a consequence, the spectral currents are given by linear circuit theory 
rules where local voltages are replaced by the electron distribution functions 
/ at the nodesB^I The spectral conductance g{e) is the average conductance 
for electrons with a given energy e : 



9ie)=S[l —]-^ = aMS[l ^. , J -i=G^ + ^g(g) (9) 



/ —J =^nS[ — -J— 

/o cr Jo cosh d2{e,x) 

Here, we do not consider the change in conductance of the opened channels 
with increasing voltage. This was included in the framework of the scattering 
matrix theory by Lesovik et alBll 

Let us consider the generic example of a sample made of a quasi- ID N 
wire "n" of length L between a S island and a N reservoir, see Fig. |2| inset. 
The spectral conductance is a measurable quantity, since it coincides with 
the zero-temperature differential conductance at voltage V = e/e : 

9{e) = ^{V = e/e,T = 0) (10) 

The spectral conductance g{€) is strongly energy-dependent and always larger 
than the normal-state conductance. In the same way, the excess conductance 
6g{e) is always positive. At finite temperature, the differential conductance 
is the integral of the spectral conductance multiplied by the energy derivative 
of the Fermi distribution function : 

This is equivalent to averaging the spectral conductance over an energy 
window of width AkBT. 

Real transport experiments involve more complex circuits than the 
generic two-probe configuration considered above. Let us consider the gen- 
eral case of a network of ID wires connected to several normal reservoirs 
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and a unique superconducting electrode. The chemical potential of the su- 
perconducting island is taken as the reference potential {V = 0). One can 
then define a spectral conductance matrix g for the whole circuit connect- 
ing the spectral currents matrix i incoming from the reservoirs with the 
non-equilibrium distribution functions / :E3 

! = £•/ (12) 

Integration of this set of relation over the energy gives a direct access to the 
current-voltage characteristic, since the summation of the spectral currents 
leads to the measurable current while integration of the / functions of the 
reservoirs leads to the chemical potential of the reservoirs. Let us note that a 
zero total current does not mean that the spectral current should be zero. In 
a reservoir with no injected current, the low-energy and high-energy spectral 
currents are in general different from zero but cancel out once integrated over 
the whole energy spectrum. 



2.3. The re-entrance effect 

As an example for the calculation of the spectral conductance, let us 
come back to the generic sample geometry made of an N wire between an S 
island and an N reservoir. The numerical solution of the Usadel equations for 
the spectral conductance g{€) in units of the normal-state conductance Gn 
is shown in Fig. g. From the calculation, the spectral conductance shows a 
maximum of about 1.15 Gn at an energy close to 5.1 ec with ec = fiD / L? . At 
higher energy, the spectral conductance decays as l/\/e. The most striking 
result is that at zero energy, the normal-state conductance is recovered. This 
is the re-entrance effect. 

The occurance of the re-entrance of the metallic conductance does not 
depend on the precise 
quasiclassical theory,c3'l 



pmetry. It is an exact resiilLfrom the 
1 the random-matrix theoryllaEl and the 
Bogoliubov-de-Gennes equations.E3 The zero-energy conductance of a prox- 
imity superconductor coincides with the normal-state conductance only in 
the absence of electron-electron interactions. In the presence of interactions 
in the normal metal, the zero-temperature conductance is predicted to differ 
from the normal-state value Gm- The zero-temperature conductance should 
decrease in the case of repulsive interactions and increase in the case of 
attractive ones £3 The ferromagnetic metals provide a case where the inter- 
actions are expected to play a great role in the behaviour of the proximity 
effect .H 

Let us try to draw a physical picture for the re-entrance effect. The 
conductivity peak for electrons with a given energy e appears in the vicinity 
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Fig. 2. Energy dependence of the spectral conductance g{e) in units of the 
normal-state conductance Gn for a N-S sample. The sample is made of a 
normal metal wire "n" of length L between a N reservoir and a S island, see 
inset. The spectral conductance was calculated with the help of the full non- 
linear Usadel equations (full line) and in the linear approximation (dashed 
line). The gap and the phase-breaking length were assumed to be infinite. 
The N-S interface transparency is considered as equal to 1. 



of the distance L^ from the S interface. Strikingly, it is the point where 
the Andreev pairs break. At smaller distances x < L^, the Andreev pair 
is a closed object, i.e. the electron and the hole follow exactly the same 
path. Another conduction electron cannot enter the Andreev pair. The 
local conductivity is unchanged by the proximity of the superconductor. 
The density of states is nevertheless depleted since the probability to find 
an isolated electron is very small. At a distance x ~ L^, the electron and 
the hole have distinguishable trajectories. The Andreev state can couple the 
electron with a hole a Fermi wave-length away. The hole can even interact 
again with the superconducting condensate and be reemitted as an electron. 
This "delocalization" of the electron enhances the local conductivity. At 
high energy, the conductance enhancement decays because the electron and 
the hole are increasingly decor related. 
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2.4. The linearization 

For the easiness of presentation and understanding, it may be convenient 
to linearize the Usadel equation and consider the regime of an infinite phase- 
breaking length. This gives : 

nD5le + 2ie9 = (13) 

The linearization is fully justified only in the case of a (thin) tunnel barrier 
between N and S. In all cases, it remains a good approximation at a distance 
from the interface larger than L^. 

The linearized Usadel equation is a diffusion equation for the pair am- 
plitude in the real space dimension x. The characteristic diffusion length is 
the coherence length L^ = yJfiD je. At a given position x, the character- 
istic energy scale is ex = fiD /x^ which coincides with the Thouless energy 
ec = hD /L? in the case x = L. 

Let us now consider some ideal geometries. 

In the case of an semi-infinite length N wire in contact with a S island, 
the angle has the following simple form : 

e = |exp[(z-l)-^] (14) 

The local conductivity can be analytically written as : 

o-(e,x) =aArcosh^[-exp(-— )sin— ] (15) 

In the case of a N wire of finite length, the N reservoir imposes a zero 
value of 9 at the boundary. Fig. |2| shows the energy dependence of the spec- 
tral conductance g{e) in units of the normal-state conductance Gn- The 
curves are calculated without and with the linearization approximation. We 
observe that the spectral conductance is qualitatively the same in the two 
cases. Although the linearization does not describe accurately the contri- 
bution of low energies, the integrated conductance is very well described. 
The position and the value of the maxima of conductance are very close 
in the two calculations. This comparison is the justification for the linear 
approximation used in Ref. lj. 



2.5. Influence of the physical parameters 

Let us consider the effect of several physical parameters on the spectral 
conductance of a sample made of a "n" metallic wire between one N reservoir 
and one S island, see Fig. inset. 
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Fig. 3. Energy dependence of the spectral conductance g{e) in units of the 
normal-sate conductance Gn for various values of the phase-breaking length. 
Each curve is labeled with the related ratio of the phase-breaking length L^ 
over the sample length L. The gap of S is taken as infinite. The barrier 
transparency is assumed to be 1. 



2.5.1. The phase-breaking length 



Up to now, we considered that the phase-breaking length L^ is much 
larger than the sample length L. In this case, phase-breaking events have no 
effect on the spectral conductance since it is the N reservoir which limits the 
diffusion of pairs in the normal metal wire. Now let us consider the opposite 
case where the phase-breaking length is smaller than the sample length. 

Fig. m shows the calculated spectral conductance for the Fig. ^ sample 
geometry but with a varying value for the phase-breaking length L^. As 
soon as L^ becomes of the order of L, the spectral conductance is affected 
at energies below e^ = hD/L'^. When L^ < L, the position and absolute 
amplitude of the spectral conductance maximum do not vary with the sample 
length L anymore, but depend directly on the phase-breaking length. The 
spectral conductance maximum shifts to 1.2 e,^ with an amplitude equivalent 
to a 16 % increase for the conductance of the length L^ of the wire. 

The phase-breaking brings a very strong cut-off for the proximity effect 
so that it plays the role of an effective sample length. As a consequence, 
the energy e^ also replaces the Thouless energy €c- In a given sample, the 
phase-coherence length can be shortened by the flux induced in the width w 
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Fig. 4. Energy dependence of the spectral conductance g{e) for various 
values of the transparency of the barrier at the N-S interface, in units of the 
normal-state conductance of the N metal alone. The calculated conductance 
takes into account both the metallic conductance of the N metal wire and 
the finite conductance of the N-S interface. The curves are calculated for 
different values of the ratio Lt/L of the barrier equivalent length Lt over 
the sample length L. The gap of S is taken as infinite. The phase-breaking 
length is considered as infinite. 



of the N wires by an applied magnetic field H. The effective phase-coherence 
L^p follows the relation xB 



L-\H) = L~\Q) + 



3 $2 



(16) 



where 00 = h/2e is the flux quantum and Lip{0) is the zero- field phase- 
breaking length. At zero magnetic field, the characteristic energy scale of 
the proximity effect is the minimum of the Thouless energy €c and the phase- 
breaking-related energy e^. As the magnetic field is increased, the energy 
€^p is decreased and the temperature of the conductance maximum is shifted 
to higher temperature. This has been observed in a previous experiment on 
the re-entrance effect .113 
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2.5.2. The interface transparency 

In this study, we are mainly interested in the case of a clean, metallic, 
interface between the normal metal and the superconductor. Nevertheless, 
it is worth knowing how much an intermediate interface transparency can 
modify the spectral conductance. The transparency t of the N-S interface 
can be smaller than 1 due to the presence of a small tunnel barrier at the 
interface. A mismatch between the Fermi velocities of the two metals can 
also induce a finite reflection coefficient at the interface. 

The interface transparency can be implemented in the calculation as a 
different boundary condition for the angle 6 at the N-S interface : 

sm{es-eN) = Lt5J (17) 

where the length Lt = le/t is the barrier equivalent length,c3 O^q and 6$ are 
the values of the angle on the N or S side of the N-S interface, respectively. 
The length Lt has a simple physical explanation in the limit of a small 
transparency since it is the length of normal metal wire which has the same 
resistance as the barrier. 

Fig. ^ shows the conductance of the total N-S sample, including the 
interface resistance. As already discussed in Ref. cJ, it shows a cross-over 
between two distinct regimes as the interface transparency is decreased. At 
large transparency we observe a maximum of the spectral conductance at 
finite energy : it is the re-entrance effect. When the interface transparency is 
so small that the length Lt is larger than the sample length L, the maximum 
of spectral conductance sits a zero energy. This is the zero-bias anomaly 
of a N-I-S tunnel junction conductance. The zero-bias anomaly has been 
first observed by Kastalsky et al.li2l It can be modulated by a flux in a 
loop geometryEB The physical interpretation is that the tunneling of pairs 
through the tunnel barrie r is enhanced by the confinement by the disorder 
in the N-metal layer.Blii0 



2.5.3. The superconductor gap 

If both the gap in S and the phase-breaking length are considered as 
infinite, the boundary condition for the complex angle at the N-S interface 
is as simple as : ^ = i7r/2. This assumption is valid at low bias and in a 
restricted temperature range, i.e. not too close to the critical temperature 
of the superconductor. Taking into account the finite value of the S gap, the 
boundary condition for the angle 9 has to be changed to : 

= |+iargth(^) (18) 
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Fig. 5. Energy dependence of the spectral conductance for various values of 
the superconducting gap in S. The curves are calculated for different values 
of the ratio A/cc of the superconducting gap A over the Thouless energy ec- 
The barrier transparency is assumed to be 1. The phase-breaking length is 
considered as infinite. 



Fig. ^ shows the calculation results for various values of the ratio be- 
tween the Thouless energy and the S gap. If one considers the whole be- 
haviour of the spectral conductance g{e) as a function of the energy e, the 
infinite gap assumption is valid only in the regime ec <C A. When the gap 
becomes of the order of the Thouless energy, a peak of the spectral conduc- 
tance enhancement is observed at e ~ A in addition to the re-entrance peak 
at e ~ 5ec- This additional peak is due to the singularity of the density 
of states in the superconductor at the gap energy. The amplitude of this 
additional peak in the spectral conductance can be very large. This peak 
has been recently discussed in Ref. c3. 

At low values of the gap A < 20 ec, only one peak at e ~ A is visible. 
Let us note that even at A = 20 ec, the contribution of the peak of the 
spectral conductance near A cannot be neglected. This will be confirmed in 
the experimental discussion. 



THE EXPERIMENTAL RESULTS 



The spectral conductance and the re-entrance effect 




Fig. 6. Micrograph of the sample. The sample is made of three arms of Cu 
wire joining the two Cu reservoirs "1" and "2" and one Al superconducting 
island "3". Each arm is 200 nm long, 80 nm wide and 50 nm thick. The 
voltage probes of the two Cu reservoirs are visible. The Al island is also 
connected to two measurement probes. For clearity, a thin white line has 
been drawn around of the Al island. Transport measurement were carried 
out both between the two Cu reservoirs ("N-n-N" geometry) and between 
one Cu reservoir and the Al island (N-n-S geometry). 



3.1. The sample configuration 

Fig. y shows the micrograph of the sample we designed for measuring 
the re-entrance effect. Preliminary results from this sample were previously 
published in Ref.123. This T-shaped sample is made of three Cu arms joining 
three wide banks : two are made of Cu, one is of Al. The Al island is 
superconducting below about Tc = 1.3 K. The Cu-Al interface has been 
carefully prepared in view of obtaining a highly transparent interface. From 
the fit to the theory (see below), we derived a maximum resistance value 
of 2 for the N-S interface of area about 200X80 nm^. This corresponds 
to a length Lt of about 80n?7T-, which is of the order of the length of the 
overlapping region between the Cu and Al layers. 

The conductance can be measured between the two Cu reservoirs ("N- 
n-N") or between one Cu reservoir and the Al island ("N-n-S"). With the 
first geometry, the redistribution of current paths in the vicinity of the N-S 
interfaced due to the superconductivity of Al has a reduced effect on the 
measured conductance. The normal-state conductance measured between 
the two Cu reservoirs is Gi2,Ar = 0.099 5". This gives a mean free path l^ 
of Gnm and a diffusion coefficient D = vple/S = 30 crra^/s, using pl^ = 
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1.5710^771.5^^ for Cu. From this data, we can calculate a value for the 
Thouless energy €c = 12/iey associated with the Cu wire length between 
the two Cu reservoirs. This value is much smaller than the gap A ~ 190 ^eV 
of Al. 

The normal-state conductances of the left (labeled 1) and right (2) arm 
were measured to be Gin = 0.190 5" and G2N = 0.206 5 respectively. The 
conductance of the upper arm in direct contact with the Al island could not 
be measured directly since the voltage probe was not in close vicinity of the 
N-S interface. From the fit of the experimental results to the theory, we infer 
a value G^n = 0.183 5 which is compatible with the sample geometry. In 
the following, all the transport measurements were carried out with an ac 
bias current corresponding to a voltage modulation of 2ij,eV. This results 
in a temperature smearing of 23 mK. 



3.2. The transport measurements 

Fig. M shows the temperature dependence of the "N-n-N" conductance 
Gi2- As the temperature decreases below the superconducting transition 
of Al, the conductance first rises. It afterwards reaches a maximum near 
T = 450 mK and eventually decreases down to the lowest temperature 
(40 mK). This behaviour is characteristic of the re-entrance effect. Here, the 
conductance decrease amplitude is comparable to the amplitude of the con- 
ductance increase at higher energy. Nevertheless, the low-temperature limit 
of the conductance is significantly higher than the normal-state conductance. 

The "N-n-S" conductance G13 measured between the Cu reservoir 1 and 
the Al island is shown in Fig. ^. The behaviour is qualitatively similar but 
with quantitative differences. The conductance enhancement is about 10% 
between 1.2 K and 40m,K, instead of 2.5% for G12. In the vicinity of the 
critical temperature of Al, the conductance drops sharply. This behaviour is 
difficult to analyze because the voltage drop in Al measured in series is also 
expected to depend strongly on the temperature. 

The bias dependence of the non-linear "N-n-N" conductance was mea- 
sured at 40 mK with an ac modulation of the bias current superposed to a 
dc current bias. Data shown in Fig. ^ exhibit a behaviour which is very 
similar to the temperature dependence data, both in energy dependence and 
amplitude. The differential conductance is minimum at zero bias voltage. It 
shows a maximum at a finite bias voltage of about 40 fieV and eventually a 
decrease at large bias. 

In summary, the conductance of N-S structure is shown to be maximum 
at finite temperature or bias voltage. These two features bring the proof for 
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Fig. 7. Temperature dependence of the "N-n-N" sample conductance G12 
measured between the two N reservoirs at zero bias. The experimental curve 
(Exp.) is shown in parallel with two theoretical curves (Th.) calculated with 
the non-linear Usadel equations. The dotted line has been calculated with 
the assumption of an infinite gap in S and the nominal ec = 12 ^eV. The full 
line is with Cc = 15.5 fieV and takes into account a constant value for the gap 
in the Al superconducting island. The phase-breaking length is considered 
as infinite and the interface transparency is taken as equal to 1. 



the re-entrance of the spectral conductance since both experiments probe the 
spectral conductance. In the zero bias and variable temperature experiment, 
energy is driven by the bath temperature (e ~ /cbT). In the very-low- 
temperature and variable bias experiment, energy is driven by the chemical 
potential of the reservoirs. 

The re-entrance of the metallic conductance in a mesoscopic proximitv 
superconductor has been previously observed in other metallic samplegi^llZI 
and in semiconductor structures.t3 The reentrance of the magnetoresistance 
oscillations has also been tracked as a function of the bias voltage in semiconductor- 
superconductor structures, the semiconductor being a two-dimensional elec- 
tron gas.EJ 



DISCUSSION 
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Fig. 8. Temperature dependence of the "N-n-S" sample conductance G13 
measured between one Cu reservoir and the Al island. Same experimental 
conditions and fit parameters than previous Figure. 



4.1. The detailed spectral conductances 

We used the quasiclassical theory to quantitatively describe the experi- 
mental results. First, the non-linear Usadel equation (Eq. ^) was numerically 
solved in order to obtain the function 92{x, e). Here we used ideal boundary 
conditions : 9 = at the N reservoir, perfect N-S interface with a step-like 
energy gap (no self-consistency), continuity for 9 and 5x9 at the central node. 
The three spectral conductances gi{€), (72(e) and 53(e) are afterwards derived 
from the ^2 function. 

Fig. ^ shows the three theoretical spectral conductances of the three 
Cu arms of the sample. The energy unit is the Thouless energy Ec = hD/L^ 
with the length of each branch assumed to be equal to L/2. The Thou- 
less energy was taken as a fit parameter. In the following, we will see that 
the agreement between theory and experiment is better with a Thouless 
energy of 15.5fj,eV. This value is indeed compatible with our experimen- 
tal uncertainty on the physical dimensions of the sample. The gap of the 
superconductor Al is A ~ 190fieV. 

From Fig. ^, the spectral conductances of wires 1 and 2 are maximum 
at 3.8 ec- The amplitude of the conductance enhancement (here 4.8%) is 
significantly smaller than in Fig. ^ because the S island is separated from 
the wires 1 and 2. A slight local maximum at higher energy is related to 
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Fig. 9. Voltage dependence of the differential conductance of the same sam- 
ple in the N-n-N measurement geometry. The theoretical curve labeled (QC 
Th.) is the result of the calculation assuming perfect Cu reservoirs, i.e. with 
a electron Fermi distribution at the phonon bath temperature. The curve 
labeled (QC Th. + res.) includes the additional effect of heating of the 
reservoirs by the bias current. The fit parameters are the same as in the two 
previous figures. 

the gap of Al. The conductance of the arm 3 is much more affected by the 
finite amplitude of Al gap. This is because Andreev pairs with an energy 
close to the gap remain coherent only very close to the N-S interface. The 
maximum spectral conductance of arm 3 is found slightly below the gap of 
Al at energy e = llec ^ A. The amplitude of this peak {Sg^/GsN ~ 100%) 
is very large. This is remarkable, since the gap of S only appears in the 
boundary condition at the N-S interface. 



4.2. Prom the spectral conductances to the I-V characteristics 

Following Eq. |l^, a matrix relation makes the link between the spectral 
current in each of the three wires and the out-of-equilibrium function / at 
the nodes of the sample. In our 3-branches circuit, we have : 



h = ai [h - In] 
h = 92 [/2 - In] 
h = -93 In 



(19) 
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Fig. 10. Detail of the calculated spectral conductances of each of three arms 
constituing the T-shaped sample. The assumed normal-state conductances 
are Gin = 0.190 5, G2N = 0.206 5 and G^n = 0.183 5. The interface 
transparency is equal to 1 and the phase-breaking length is considered as 
infinite. The gap A of S is equal to 12 times the Thouless energy. 



where all the quantities are a function of e. Here /i, /2 and /at are the 
out-of-equilibrium part / of the energy distribution functions in reservoirs 
1, 2 and at the central node. Since the node 3 is the superconducting island, 
its out-of-equilibrium distribution function /s is zero. The superconducting 
island 3 is the reference for the voltage. In a N- metal reservoir at voltage V 
at perfect thermal equilibrium at temperature T, f is of the form EB 






(20) 



According to the Kirchoff law ii + ^2 + is = at each energy and following 
Eq. |l^, /tv is a linear combination of the reservoirs distribution functions : 

51 /i + 92 /2 



/ 



N 



91+92+ 93 
where again all the quantities depend on the energy e. 
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The I-V characteristics were calculated by integration of the equations 
19| taken into account the boundary conditions. In the two experimental 

(22) 



configurations, they are : 

N-n-N : h = -/, h = I, h = 0, V = Vi - V2 



N-n-S : h = -I, h = 0, h = -/, V = Vi 

Let us first consider the zero-bias temperature dependence of the con- 
ductances. In both Fig. and ^, calculated curves are shown in comparison 
of the experimental one. Calculated curves are shown in the two cases of an 
infinite gap taking ec = 12 iieV and with a finite gap A = 12ec = 186 /ieV 
while Ec = 15.5 fieV. In the last case, the Thouless energy corresponds to a 
Thouless temperature ec/ks = 180 mK. The ratio A/ec is about 12. 

The agreement between experiment and theory is satisfactory in the 
N-n-N geometry in both cases. In the N-n-S geometry, the infinite gap 
assumption does not provide a good description of the data. The agreement 
is good in the case of a finite gap and ec = 15.5 ^ueF. The fit is not very 
accurate near the critical temperature Tc. This is because our assumption 
of a constant value for the superconducting gap is not valid anymore in this 
region. 

From these fits, we conclude that the experimental temperature de- 
pendence data are well described by the theory. The value of the physical 
parameters introduced in the calculation are within the measurement uncer- 
tainty. The effect of the finite gap is non-negligeable, especially in the N-n-S 
geometry where the conductance of arm 3 contributes. This is because the 
pair amplitude F(e, x) is very large at energies close to the superconductor 
gap. The Andreev pairs with this value of energy remain coherent only very 
close to the N-S interface. 

At zero energy and zero temperature, one should get the normal-state 
value of the conductance. This is not seen in th^^experiment. This may be 
an effect of the interactions in thenormal metal.c3 However, our experiment 
cannot bring a definitive answer .lHI 

Let us now consider the finite bias data. As expected the re-entrance 
effect, i.e. a conductance peak at energy close to the Thouless energy, is 
observed in both the calculated and measured data. However, there is no 
quantitative agreement, as the peak position is at significantly lower energy 
than expected. 



4.3. The role of electron energy distribution in the reservoirs 

In order to obtain a better description of the bias dependence experi- 
ment, one has to consider the broadening of the energy distribution function 
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in the Cu reservoirs by the injected current. We choose to describe this heat- 
ing effect by an effective temperature T(,ff of the reservoir. This effective 
temperature is different from the phonon temperature and depends on the 
bias current which is injected in the reservoir. 

Two different physical effects may be the cause for a heating of the reser- 
voirs. First, the chemical potential may be ill-defined in the reservoirs. If the 
reservoir resistance per unit length p^ is not sufficiently small, a significant 
potential drop piLinl appears in the reservoir, where Lin is the inelastic 
mean free path. The electron energy distribution will be close to a Fermi 
distribution with an effective temperature T^jf. If the inelastic mean free 
path is temperature- independent, then the effective temperature behaves as 
Teff oc /. 

A second possibility is the Joule effect induced in the reservoirs due to 
the bias current passing through. This heating power will be evacuated by 
the phonons. The phonon temperature T will be again different from the 
effective electron temperature Teff. Let us first consider the dissipation due 
to local Joule heating throughout the sample. The effective teinDerature 
then expresses as : T^fj ex /^'^ in the limit where T^fj » T.EJ In the 
pure mesoscopic case where the dissipation occurs only in the reservoirs, the 
Joule effect due to the sample resistance will be evacuated in the reservoirs. 
The reservoir temperature becomes therefore non-uniform and the diffusion 
equation for heat has to be solved. Again in the limit where T(,ff >> T, the 
exponent of the previous power-law will be changed into 4/7. 

The current modulation of the sample creates a modulation of the ef- 
fective temperature T^ff. Then, the derivative of the electron energy distri- 
bution function writes : 

^=4 + ^'^ (23) 

dV de dieff dV 

This introduces an additional term to the usual integration of the spectral 
conductance over a thermal window, namely a convolution of the spectral 
conductance with the derivative of a peak function. As a consequence, the 
usual term of the conductance peak is shifted in the experiment compared 
to the case of ideal reservoirs. 

We described the energy distribution of the electrons in the reservoirs by 
a Fermi distribution function with a current-dependent effective temperature 
Tgjj. The effective temperature follows a power law : T(,ff = CI", with C 
being a constant common for all the data. We found a better agreement by 
choosing a = 1 as compared to the other possible exponents. The ajusted 
value of the constant C gives an estimate for the inelastic mean free path 
Lin — ^pm. The result of the calculation is shown in parallel with the 
experimental curve in Fig. H. The agreement between experiment and theory 
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is quite improved. The position and amplitude of the conductance peak are 
well described. 



4.4. A direct test of the non-equilibrium distribution function 

Here we describe a complementary experiment designed to test the as- 
sumption of independent energy channels. The strong point is that along 
the wire, the electron distribution function is strongly non-equilibrium, see 



e.g. Eq. 21. Indeed, it consists of a superposition of step functions centered 
at the chemical potential of the respective normal reservoirs. The electron 
energy distribution has been directly measured in Ref. e3. 

The concept is to use the part 3 of the sample as a probe for the dis- 
tribution function at the node of the sample. The part 3 of the sample has 
a strongly energy-dependent conductance, like a N-I-S tunnel junction. We 
injected a dc current /„ from the Al probe (3) to the left Cu probe (1) while 
we measured the resistance of the left arm (1) of Cu. This was made by 
biasing with a low- frequency ac current between the two Cu probes (1 and 
2) and measuring the voltage drop between the Al probe (3) and the left Cu 



probe (1), see Fig. 11a inset. With a zero injection current, the data show 
the same behaviour as the full Cu wire, see Fig. ^. With a non-zero injection 
current, the differential conductance becomes asymmetric in respect of the 
bias current. Peaks arise, which shift with varying the injection current from 
4 to 10 fiA. 

In this configuration, the current in the three branches are /i = — (/ + 
Iq), I2 = I, I3 = Iq- As already stressed, the part 3 has the strongest en- 
ergy dependence of the spectral conductance. For simplicity, our qualitative 
discussion will neglect the proximity-induced enhancement of the conduc- 
tance of branches 1 and 2 and keep the contribution of the branch 3 only. 
We will restrict our qualitative discussion to the zero-temperature limit. By 
integrating and inverting Eq. ^, we obtain : 

yiN LrlAT <~r3N 

V2 = -^-^ + 6Vn (24) 

The non-linear contribution SVn contains the effect of the energy-dependent 
conductance enhancement Sg^^e) : 

SVn = - r ^-§^fN{e)de (25) 
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Fig. 11. (a) : Measured differential conductance of the left arm of the 
sample as a function of the bias current for various values of the hot electron 
injection current Iq. Temperature is 100 mi^. (b) : Calculated differential 
conductance at T = 100 rriK of the left arm of the sample as a function of 
the bias current for various values of the injection current Iq. In both curves 
sets, the dashed line indicates the normal-state conductance. 
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In our experiment, we measure the differential resistance dVi/dl versus 
/ with /q kept constant by an independent floating source. In the normal- 
state 1/GiN is measured. Below Tc, the non-linear term brings a significant 
contribution. It measures directly the important features of the distribution 
function at the central node. 

dVi. 1 r ^9z{e)dfN{e) 



dl CiTv J~oo G^N dl 

At very low temperature, the derivative of the distribution function of 
the reservoirs Vi and V2 should tend to 5— functions centered at the chemical 
potentials of the reservoirs. In this hypothesis, it is straightforward to trans- 
form the convolution product into a combination of Sg3{eVi) and 5g3{eV2) 



using Eq. ^, 19 and 24. The final formula for the conductance dl/dV (here 
V = -Vi) is : 

^1/0 = G^^i^ + ^^i^93{eV2) - dgsieV,)}] (27) 

with G]\f = Gij\f + G2N + G37V- 



From Eq. ^, we notice that both Vi and V2 are tuned by the in- 
dependent driving currents / and Iq. They vanish respectively for I~ = 
— /o(l -t- Gin/G^n) and 1+ = /q G^n /G2n- Referring to Fig. |l^, the strong 
minimum of (5(73 (e) at e = show up as a minimum of ^^ at I^ where 
V2 = and a maximum at /~ where Vi = 0. This is precisely the origin of 
the strong assymetric behaviour observed in the experimental curve of Fig. 



ll| . The asymmetry is present because of the distribution function J^n at 
the central node exhibits sharp steps at the chemical potentials of the reser- 
voirs eVi and eV2- With this experiment, we have proven that the energy 
distribution of the electrons at the node in the middle of our sample is truly 
out-of-equilibrium. In brief, we are indeed in a mesoscopic regime. 

In order to achieve a quantitative comparison between theory and ex- 



periment, we solved the full matrix expression Eq. 12 which makes no simpli- 
fication on 5gi{e) and 5g2{e).We included the heating effect in the reservoirs 
in the same way than in the re-entrant conductance data. The data of Fig. 
|ll| are well described by using the same constant C and exponent 1 as in the 
Fig. 1^ and ^. As expected, the two extrema occur at values of I proportional 
to Iq. When Iq = 0, this contribution is absent. Only the smaller symmet- 
ric contribution of dgi (e) remains. Our previous analysis leading to Eq. ^ 
ignores this contribution but highlights the main contribution due to 5g3{e). 
We als(i_studied an alternate geometry with the inversion of the currents / 
and /q-" In this case, the same two terms Sgs{eVi) and Sg3{eV2) are added 
instead of being subtracted. The experimental data (not shown) were again 
consistently fitted by the calculation with the same physical parameters. 
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5. CONCLUSION 

In conclusion, we presented a thorough study of the spectral conduc- 
tance of a mesoscopic normal metal in contact with a superconductor. The 
spectral conductance exhibits a re-entrance effect at zero energy, i.e. the 
zero-energy spectral conductance coincides with the normal state conduc- 
tance. From the theoretical calculations, we observe that the spectral con- 
ductance is sensitive to the absolute value of the many physical parameters 
including : phase-breaking length, gap of the superconductor and interface 
transparency. 

We performed experimental measurement of the reentrance effect as a 
function of both bias voltage and temperature. The experimental data is 
well described by the quasiclassical theory. Nevertheless, the description of 
the non-linear conductance data requires taking into account of the heatinff 
of the N reservoirs by the bias current. Compared to a previous study,li3 
we were able to describe quantitatively the reentrant resistance with the 
quasiclassical theory. Thanks to the well-controlled geometry of the normal- 
metal reservoirs, no scaling factor was necessary. 

The large energy-sensitivity of the conductance enhancement by the 
proximity effect can be used to probe directly the energy distribution func- 
tion without using a tunnel junction. Thus we were able to measure the con- 
ductance enhancement and test the distribution function in a single sample. 
The whole set of data agree with the theoretically-calculated curves with 
a single set of physical parameters. This confirms that most of the be- 
haviour of this sample is understood by the theory of the proximity effect in 
a non-interacting metal. In contrast, the difference between the conductance 
measured in the zero-temperature limit and the normal-state value raises the 
question of the importance of the effect of interactions. 

We thank F. W. J. Hekking, T. Martin, B. Spivak, A. F. Volkov and F. 
Zhou for fruitful discussions. This work was supported by the DRET, the 
Region Rhone- Alpes and the TMR contract "Dynamics of superconducting 
nanocircuits" from EU. 
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